Scalar and vector decomposition of the nucleon self-energy in the 

relativistic Brueckner approach 

C. Fuchs, T. Waindzoch, Amand Faessler, and D.S. Kosov 

Institut fiir Theoretische Physik der Universitdt Tubingen, 

D-72076 Tubingen, Germany 
We investigate the momentum dependence of the nucleon self-energy in 

nuclear matter. We apply the relativistic Brueckner-Hartree-Fock approach 
and adopt the Bonn A potential. A strong momentum dependence of the 
scalar and vector self-energy components can be observed when a commonly 
used pseudo-vector choice for the covariant representation of the T-matrix 
is applied. This momentum dependence is dominated by the pion exchange. 
We discuss the problems of this choice and its relations to on-shell ambigui- 
ties of the T-matrix representation. Starting from a complete pseudo-vector 
representation of the T-matrix, which reproduces correctly the pseudo-vector 
pion-exchange contributions at the Hartree-Fock level, we observe a much 
weaker momentum dependence of the self-energy. This fixes the range of the 
inherent uncertainty in the determination of the scalar and vector self-energy 
components. Comparing to other work, we find that extracting the self-energy 
components by a fit to the single particle potential leads to even more am- 
biguous results. 
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I. INTRODUCTION 

Concerning the nuclear many-body problem the relativistic Dirac-Brueckner-Hartree- 
Fock (DBHF) approach turned out to be remarkably successful in describing the nuclear 



matter saturation mechanism. The major improvement compared to non-relativistic treat- 
ments is based on an additional density dependence introduced in the formahsm by using a 
self-consistent spinor basis. To solve the Bethe-Salpeter equation, i.e. its three-dimensional 
reductions, a variety of approaches relying on different techniques and various bare nucleon- 
nucleon interactions have been developed over the last decade P|-pi]|. All these calculations 
are able to describe reasonably well, although not excellent, nuclear matter properties. In 
the meantime two main results can be regarded as settled: first, the nuclear single parti- 
cle potential originates from the cancelation of large repulsive vector and attractive scalar 
fields, and second, the magnitude of the effective mass is reduced to a value of about ~ 0.6M 



at saturation density. Both findings are consistent with effective hadron field theories |11 



where, e.g. the effective mass can be determined from the spin-orbit interaction in finite 



nuclei 121. Refined treatments which take into account hole-hole excitations 1131 or include 



an additional scaling of meson properties in the nuclear medium [|T^ do not significantly 
alter these results. 

To test the DBHF approach over a wider range of physical problems also finite nuclei 
have been studied within effective Hartree-Fock calculations |T5| and within a density de- 
pendent hadron field theory P^JT7| . For a successful apphcation of DBHF results to heavy 
ion collisions using relativistic transport models, however, it is required to go beyond a local 
density approximation and to properly account for the non-equilibrium features of the highly 
anisotropic phase space configurations in such reactions |T^. As soon as nuclei overlap and 
particles are positioned at close quarters in configuration space, the momentum dependence 
of the nuclear self-energy starts to play a crucial role in the description of heavy ion collisions. 
Unfortunately, the determination of the momentum dependence of the nucleon self-energy is 
a subtle problem in the DBHF approach which has not yet led to settled results. Generally, 
the techniques applied in the standard relativistic Brueckner approach rely on a weak mo- 
mentum dependence of the self-energy inside the Fermi sea. This assumption was supported 
by various calculations in the past |T|-§[. In Ref. 0] it was noted that the determination of 
the self-energy leads to ambiguities arising from the projection of the T-matrix onto positive 



energy states and that the actual momentum dependence strongly depends on the choice of 
the used nucleon-nucleon interaction 0. These ambiguities can be avoided when negative 
energy states are taken into account in the calculation However, the conventional 

nucleon-nucleon potentials are determined for particle-particle scattering and an extrap- 



olation to anti-particles is in itself ambiguous. To avoid the latter ambiguity we will restrict 
our discussion merely to the particle sector. 

In practice there exist various procedures to determine the self-energy. In the method 
proposed by Horowitz and Serot [|l|] one projects onto Lorentz-invariant amplitudes of the 
T-matrix (or in-medium G-matrix) This method is, however, not unique since 

pseudo-scalar (PS) and pseudo-vector (PV) matrix elements are equivalent for on-shell par- 
ticles in the positive energy sector. Hence, it is impossible to disentangle pseudo-scalar and 
pseudo-vector contributions in the on-shell T-matrix. On the other hand, the pion exchange 
contribution is known to react rather sensitively on the particular choice of its representation 
pl| , p!l|] , i.e. PS or PV. Thus in the past different choices for the representation of the cor- 
responding amplitude have been used in the relativistic Brueckner scheme [00,0. We find 
that these choices strongly influence the magnitude and in particular the momentum depen- 
dence of the self-energy in the nuclear medium. In the present work we will discuss various 
possibilities to decompose the T-matrix which are based on the general representation of 
covariant amplitudes proposed in Refs. and . A strong momentum dependence of the 
self-energy is found to originate from pseudo-scalar admixtures due to pion-exchange contri- 
butions. These pseudo-scalar admixtures are still present when the so called 'pseudo-vector 
choice' is applied, as it was done in Refs. 0,1,^0. Suppressing the undesirable admixtures 



by making use of a complete pseudo- vector representation, as proposed in Ref. [^, we find 
a much weaker momentum dependence of the self-energy components. 

In an altogether different approach pursued in Refs. [§|,|10|! the self-energy components 
are determined indirectly by a fit to the single particle potential. Here one circumvents the 
ambiguities of the projection methods. The problem is, however, only shifted to another 
level of ambiguity, since one has to extract two functions out of one. Thus, the fit method 



leads to highly ambiguous results concerning the momentum dependence of the self-energy 
components. 

The present paper is organized as follows: in section 11 we review the structure of the 
nucleon self-energy in nuclear matter. Thereafter we describe in section 111 the representation 
of the T-matrix by Lorentz invariant amplitudes and discuss the different on-shell equivalent 
choices used in the literature. In section IV we present new results obtained for the different 
choices of the T-matrix representation utilizing the Bonn A potential as the bare nucleon- 
nucleon potential. There we also discuss, as a side remark, the approach pursued in Ref. 
PJTOll. The paper ends with a summary and the conclusions of our work. 



II. NUCLEON SELF-ENERGY IN NUCLEAR MATTER 

The properties of dressed nucleons in nuclear matter are expressed by the self-energy 
which enters the in-medium nucleon propagator as the formal solution of the Dyson equation 

Due to translational and rotational invariance, parity conservation and time reversal invari- 
ance the self-energy in isospin saturated nuclear matter has the general form S = — 7^S^. 
It depends on the Lorentz invariants k"^, k-j and j^, with and A;^ being the baryon current 
and the nucleon four-momentum, respectively ||rT|. The invariants can also be expressed in 



terms of k^, |k| and kp, where kp denotes the Fermi momentum. Furthermore the vector 
part of the self energy has contributions proportional to kf^ and to the current j^. Defining 
the streaming velocity as = / the momentum k^ can be decomposed into contri- 
butions parallel and perpendicular to the streaming velocity, i.e. k^ = {k ■ u)u^ + /\^^ky 
with the projector A'^^ = g^'' — u^u'^ . The vector part of the self-energy can then be written 
covariantly as P,p^ 



S'^ = Som'^ + S.A'^'^fc^ . (2) 

Thus the full self-energy reads 



= T^sik, kp) - 7o T.o{k, kp) -k T,^{k, kp) |rf (4) 

where the subscript RF indicates the respective expressions in the nuclear matter rest frame 
{u^ = 6^^^) The Ss, So and components are Lorentz scalar functions which actually 

depend on A;o)|k| and kp. They follow from the self-energy matrix by taking the respective 
traces p2i 



Ss = \tr [S] (5) 
So = ^tr [^,u^n = ^tr [70 (6) 
Sv-^^t.[A-7.^.S]^^t.[7-kS],, . (7) 

The Dirac equation for the in-medium spinor basis can be deduced from the Green 
function ([l|). Written in terms of effective masses and momenta 

M* = M + ReJ:, , kl = kf, + Re (8) 

the Dirac equation has the form 

\^*- M* -iImJ:]u{k) = 0. (9) 

In the following we will work in the quasi-particle approximation and therefore we neglect the 
imaginary part of the self-energy from now on. Thus the effective nucleon four-momentum 
will be on mass shell even above the Fermi surface, fulfilling the relation = M*^. Since 

we only deal with the real part of the self-energy in the quasi-particle approximation we omit 
this in the notation. In the nuclear matter rest frame the four-momentum follows from Eq. 



k* = k(l + E,) , fc* = E* = ^k2(l + E,)2 + M*2 (10) 
which allows one to eliminate the S^-term in the Dirac equation, 

[a-k) --f^M*]u{k) = E*u{k) , (11) 



by a rescaling of effective mass and energy 

M* ~ E* I ~ — 

M* = , E* = = Vk^ + M*2 . (12) 

From the Dirac equation in the form (^) one derives the in-medium relativistic Hamilto- 
nian for nucleons and thereby the operator for the single-particle potential within the DBHF 
approximation, i.e. U = 7°S. The expectation value of U, i.e. sandwiching U between the 
effective spinor basis, yields the single particle potential 

which can be evaluated as 

A/f* k*Yl^ 

U(k,k,) = —Z.-^ (14) 



/k2(l + + Af*2 ^k2(l + S,)2 + M*2 

In many applications 1^,10] the single particle potential is only given in terms of a scalar 



and zero- vector component. This can be achieved in nuclear matter by introducing reduced 
fields Eg and So as 

E3 = M*-M=^^^ , E, = E*~E = E^-E*{k)E, . (16) 

Then the single particle potential has the form 

M* ~ 

[/(fc,fcp)= Ss-So . (17) 
E* 

Frequently the reduced fields, Eq. (p!6D, are used rather than the projected components 
since they represent the self-energy in a mean field or Hartree form. Thus they can easily 
be related to effective hadron mean field theory p6| , p!7| . Such a representation is meaningful 



since the Sv-contribution is a moderate correction. 

In contrast to the single particle potential which can easily be derived form the T-matrix 
1^ the extraction of the self-energy components is a subtle problem. In the latter case 
one has to represent the T-matrix within the Clifford algebra in the Dirac space which is 
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not free from severe ambiguities. Before we discuss this point in detail we briefly recall 
some basic features of the relativistic Brueckner scheme. For more details see e.g. Refs. 
|1|J^J^. The iteration of the Thompson equation requires to determine the self-consistent 
spinor basis, Eq. (^). In practice the problem is treated in terms of the reduced effective 
mass M* to avoid an explicit dependence on the space-like components ([TT]). Actually, 
the zero-vector component So does not enter the self-consistency problem. In the standard 
treatment |]l],|^,|^J|,^ the effective mass is assumed to be entirely density dependent, i.e. a 
constant effective mass M* = M*(|k| = kp,kp) generally taken at the Fermi-momentum is 
used as the iteration parameter. The effective mass at its value at the Fermi momentum 
is reinserted into the Thompson equation and this procedure is repeated until convergence 
is reached. Such a treatment is self-consistent concerning the density dependence of the 
effective interaction screened by the medium and appears to be justified if the self-energy is 
weakly momentum dependent inside the Fermi-sea. If the self-energy is, however, strongly 
momentum dependent, as it was observed in one principally has to go beyond the 

standard approach. Then one has to include the momentum dependence of the effective 
mass in the Thompson propagator as well as on the level of the interaction. For the present 
investigations we will not include such an explicit momentum dependence in the Brueckner 
scheme since first of all, this is technically a rather involved problem. Secondly, we will verify 
in the next section that the momentum dependence of the self-energy is mainly governed 
by the treatment of the pion-exchange in the T-matrix representation. Thus a more careful 
treatment of the pion-exchange leads to a much weaker momentum dependents of the self- 
energy as it is desirable for the present self-consistency scheme. 



III. DECOMPOSITION BY LORENTZ INVARIANT AMPLITUDES 



The self-energy for the nucleon k follows from the T-matrix by integrating over the 
occupied states q in the Fermi sea 



-tr 



Gri{q)T{qk; qk - kq) 



(18) 



Here T denotes the T-matrix (or G-matrix) operator depending on four Dirac indices of the 
ingoing and outgoing nucleons. Due to antisymmetrization T contains a direct ( Hartree) 
and an exchange (Fock) contribution. The Dirac propagator 

Guiq) = [f+ M*]2m5iq*^ - M*^)e{q*')e{kF - |q|) (19) 

projects onto positive energy states in the Fermi sea. In order to project out the self-energy 
components, Eqs. the T-matrix has to be decomposed into Lorentz invariants. Since 

we need to consider only on-shell scattering of particles in eq. ([I8|) , five invariant amplitudes 



with five covariants are sufficient to represent the on-shell T-matrix In this case, the 
scalar, vector, tensor, axial-vector and pseudo-scalar terms 

5=101 , V = 7^ ® 7^ , T = a^"' ®a^^ , A = 757'' ® 757^ , ^ = 75 ® 75 (20) 

form a linearly independent, however, not unique set of covariants. Using this special set, 
the on-shell T-matrix can be represented entirely by 'direct' amplitudes 

f{e) = Fi{e)s + F2{e)v + F3{e)T + F^{e)A + Fr,{d)p . (21) 

The covariant amplitudes Fi are determined from anti-symmetrized plane wave helicity 
matrix amplitudes |^ which obey the selection rule (— )^+'^+-'^ = —1. We solve the relativis- 
tic Thompson equation for the T-matrix |^ consistently in the two-particle center-of-mass 
frame. Since we need only diagonal matrix elements for the self-energy (|T8|) the Fi ampli- 
tudes are required at the scattering angle = 0. They are already anti-symmetrized and 
contain implicitly direct and exchange contributions. Therefore the simple representation 
(^) is sufficient to calculate the self-energy. If we insert the T-matrix given by Eq. (|2lD into 



Eq. (|18D and take the trace over the Dirac propagator only the scalar and vector amplitudes 
Fi and F2 survive and contribute to the self-energy, 

= - / Ki^i(O) + fFM] . (22) 

The representation of T is, however, not uniquely defined in the on-shell case and therefore 
various alternative possibilities exist to construct the set of five independent covariants in 
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the subspace of positive energies. Although the different representations discussed below 
are all equivalent if one works with the pseudo-scalar covariant P, their difference becomes 
crucial as soon as one switches from the pseudo-scalar to the pseudo-vector representation. 
The PV covariant in the medium is defined as 

PV = fcit„«fciH. (23) 

2M* 2M* ^ ' 

with k\, q\ the initial and k^., q2 the final momenta of the scattering particles. The ambiguity 
of the decomposition procedure arises from the fact that the PS and PV matrix elements 
are identical (using the Dirac eq.) for on-shell scattering of positive energy states, i.e. 

uih) iMki) = uihhMh) . (24) 

Thus the corresponding amplitudes are identical as well and it is impossible to uniquely 
disentangle the PS and PV contributions in the T-matrix. 

However, it is known from meson theory of the nuclear interaction that a pseudo-vector 
representation of the ttN coupling is preferable The influence of the pion in the rela- 



tivistic theory has been discussed in detail, e.g. in Refs. [0,^. There it was shown that the 
one-pion exchange contribution to the nuclear optical potential tends to increase drastically 
at low momenta if the vrA^ vertex is treated as PS. One reason for this behavior is a strong 
PS coupling to negative energy states which is not apparent in non-relativistic approaches. 
The PV vertex suppresses the coupling to antiparticles since the overlap matrix elements 
vanish for on-shell scattering 

Hh) {^^^ iMh) = . (25) 

Therefore a PV vertex is more consistent with the approximation scheme of the conventional 
Brueckner model where the coupling to negative energy states is not taken into account. The 
PV vertex also strongly suppresses the pion contribution in particular at low energies which 
is more in accordance with the empirical knowledge from proton-nucleus scattering ||TT],pT 



Consequently, the nN vertex of the bare interaction is usually treated as PV [20 



Due to these facts the usage of a PV covariant in the decomposition of the T-matrix 
is considered as preferable However, simply replacing P by PV in Eq. (|21| ) leads 

to identical results; firstly, because the corresponding amplitudes are equal and secondly, 
because only the scalar and vector amplitudes Fi and F2 contribute to the self-energy (p2D. 
Motivated by this fact alternative representation of the T-matrix have been used 
which lead to an - in principle superfluous - explicit splitting into 'direct' and 'exchange' 
contributions and give different results when the PS 1 — > PV replacement is performed. 



A. Symmetrized representations 
Based on Ref. p3| Tjon and Wallace discussed a representation which accounts for the 



structure of the exchange contributions in T in a more transparent way |^ . To express ex- 
change contributions also the Dirac indices of the covariants are interchanged. The transfor- 
mation S interchanges the Dirac indices of particles 1 and 2, i.e. 5*^(1)0-^(2),- = u{1)t-u{2)„. 
Thus one obtains the interchanged covariants as S* = SS, V = SV, T = ST, A = SA and 
P = SP. The interchanged covariants are related to the original covariants (pT]) by a Fierz 
transformation JF [pil] 
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(26) 



For a definite scattering angle 9 the T-matrix is now represented by five symmetrized co- 
variants 



T{e) = Me){s -s) + h{e)-{T + t) - u{e){A -a) + m{v + v) + m{p - P) . 

(27) 

These symmetrized covariants are constructed so that the amplitudes /, account for the Pauli 
principle in a transparent way. Interchanging the outgoing particles, anti-symmetrization 



10 



requires the following symmetry 



//(vr -0) = i-Y^'fm 



(28) 



for definite isospin / = 0, 1. 

The relation between the five new amplitudes /j and the former amplitudes Fi is given 



by the transformation |21 
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(29) 



With relation ( p8D one can express T as the combination of two terms which resemble a 
direct and an exchange contribution, i.e. T = — T'^, by 



f °(e) := /i(^^)5 + - /3(^^)A + U9)V + h{e)P 



(30) 



f ^(^) := (- 



A(7r-^)^ + /2(7r 



1^ 



T - /3(7r - ^^)A + U{ti - ^^)\/ + /5(7r - e)P 



However, it can not be claimed that and are the real direct and exchange matrix 
elements but only their combination yields the fully anti-symmetrized matrix elements. 
The representation proposed by Horowitz and Serot ^ and also used in other works 
is of a similar structure as the one above with the difference that it is based on direct 
amplitudes 



f^{e) := ^ [F^{9)S + F2{9)V + Fs{9)T + F,{9)A + F,{9)P] 



(31) 



^^{9) := i-Y^^l [Fi(7r - 9)S + F^in - 9)V + F3(7r - 9)f + F^{n - 9)A + F^{-k - 9)P 



To be more precise, in this approach both set of amplitudes -Fi(vr) and Fiiv: — 9) are ob- 
tained from the helicity matrix elements applying representation ( pTD at angle 9 and n — 9, 
respectively. Doing this, the amplitudes are essentially different from the fi amplitudes 
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of Eq. p7|). Instead of relation ( pSD anti-symmetrization requires now that the exchange 
amphtudes are connected to the direct amphtudes by the Fierz transformation (pG]) 



Fl{n-9) = {-Y^,.Fj{9) . (32) 



The difference between the representations ( PT| ) and (^Tj) hes in the fact that Eq. ( PT] ) 
leads to an explicit splitting into direct and exchange contributions which, however, can be 
regarded as - at least partially - artificial. Since the amplitudes Fi{9) and Fi^ir — 9) in (|3l|) 
are determined from already anti-symmetrized helicity matrix elements one has in that case 
the identity 

f ^ = -f ° . (33) 

Hence the normalization factor 1/2 which determines the splitting into 'direct' and 'ex- 
change' parts in (|3lD is judicious and could also be fixed differently by any normalized 
linear combination. In this context we want to mention that concerning the original work of 
Horowitz and Serot [|I| this statement would not hold because they used non-antisymmetrized 
(unphysical) helicity states in their formalism and therefore they had to anti-symmetrize 
explicitly the T-matrix elements by splitting the representation into direct and exchange 
contributions as done in Eq. (|3TD . As a consequence, their amplitudes Fi at angle vr and 
71 — 9 did not fulfill the anti-symmetry relation (|32D and thus only the representation (|3TD 
was physically meaningful. 

Working with physical helicity states, however, one retains relation ( ^2]) and therefore all 
types of decompositions, Eqs. (pTl-pl|), are equivalent as long as one restricts to a pseudo- 
scalar representation. But if one replaces the PS by the PV covariant the equivalence of the 
expressions (pTD, ( |27D and ( |3lD is destroyed. Our 'optimal choice' of the representation for 
using the PV covariant will be discussed in the next section. 
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B. Pseudo-vector representation 



The nucleon-nucleus potential is very sensitive to the treatment of the pseudo-scalar or 
pseudo-vector treatment of the pion-nucleon interaction. As already mentioned, the pion 
should be preferentially treated with a pseudo-vector coupling in the bare NN interaction. 
The corresponding coupling strength is fixed in such a way that it reproduces the on-shell 
PS coupling strength in the vacuum This already leads to a suppression of the vertex 
by a factor inside the nuclear medium [TT|. However, the major suppression of the 

pion exchange contribution to the nucleon self-energy originates from the different cofactors 
which arise if one inserts the Dirac propagator (|19D and the T-matrix from (^T]) or (^) into 
the Eq. ([18|) for the self-energy i.e. 

tr [{f+ M*)P] = -{f+ M*) (34) 
tr [if+ M*)PV] = (r+ M*) - . (35) 

The influence of the pion, in particular the sensitivity on the PS or PV representation, is 
most clearly demonstrated at the Hartree-Fock level. This means that T is replaced by 
the bare NN interaction V and no further medium effects are taken into account, i.e. bare 
nucleon masses are used and the Pauli operator in the Thompson equation is neglected. 
Furthermore, the Hartree-Fock expressions are known analytically which provides also 
a straightforward check of the involved projection techniques. It is a well known fact |TT 



that the PS treatment yields extremely large pion contributions to the self-energy whereas 
the PV representation suppresses these terms almost completely. This effect is demonstrated 
in Fig. 1 where the Hartree-Fock contributions of the pion only to the nuclear self-energy 
are shown. The calculations are performed for a nuclear matter density of p = 0.166 fm~^ 
with a PS and PV (denoted as 'full PV in Fig. 1) vrNN coupling of gl/4:TT = 14.9 and the 
pion form factor taken from the Bonn A interaction The bare nucleon mass is used 



in this example. It is seen that the PS description of the pion exchange yields extremely 
large self-energy components at low momenta which are rapidly decreasing with increasing 
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momentum. The Eg and So contributions are almost identical which leads to a strong 
cancelation effect in the single particle potential. The PV description (denoted as 'full PV 
in Fig. 1) suppresses the pion by nearly two orders of magnitude and even on this new 
scale the momentum dependence is much less pronounced. Remarkably, scalar and vector 
contributions have now opposite signs. Thus they add up in the potential ( [T7| ) and the 
remaining momentum dependence is not the remnant of a huge cancelation effect as in the 
PS case. 

For comparison we also show the results which are obtained in the projection scheme 
when the so-called 'PV choice' is adopted. In the 'PV choice' commonly used pjQ we simply 
replace the pseudo-scalar covariants by the pseudo- vector covariants P, P i — > PV, PV in the 
decomposition of the T-matrix (|3T|). Due to the on-shell equivalence of the respective matrix 



elements, Eq. (p^ the amplitudes remain thereby unchanged. Now the explicit choice of 
the T-matrix representation starts to play a decisive role since the PS and PV exchange 
terms contribute differently to the self-energy, Eqs. (|3^ , |35D . Thus the strength of the PS 
I — > PV replacement, determined by the respective amplitudes f^iO) and F^{tt — 6), becomes 
important. In most previous calculations the representation (^Tj) was used. The effect 
of the replacement can be seen in Fig. 1 where the Hartree-Fock contribution of the pion 
exchange to the nucleon self-energy using the PV choice is shown. First of all, it should 
be noted that the result using the PV choice for (13T|) is identical to the result which one 
obtains if one uses the replacement in the symmetrized representation, Eq. (pTf), since both 
amplitudes f5{0) and F^{7t — 6) agree in the Hartree-Fock approach with only pion exchange. 
It is obviously transparent from Fig. 1 that within the traditional PV choice the pion is 
not treated correctly as a pseudo-vector particle. The PV choice representation is rather a 
mixture of the PS and a complete PV representation. Although the pion contribution to 
the nucleon self-energy is suppressed by about a factor of two compared to the original PS 
case, the structure of the self-energies, in particular the pronounced momentum dependence 
is still very similar. The reason for this behavior of the self-energy is easily understandable. 
Even after replacing P, P by PV, PV, both representations, Eqs. (^Tf) and (^), still contain 
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pseudo-scalar admixtures because the Fierz transformation (p6D couples all covariants. This 
leads to identities for the symmetrized vector and tensor covariants |^ 

-(T + f) = S + S + P + P (36) 
V + V = S + S~P-P . (37) 

In order to completely remove the PS part from the interaction one first should use the 
identities above which leads to the following decomposition [pil] 



m = (A + /2 + h)S - (A - /2 - h)S - UA - ~A) 

+{h + f2-h)P-ih-f2 + h)P . (38) 



If the replacement P, P i — > PV, PV is now performed in Eq. (^) instead of Eq. (|27D or 
Eq. (|3TD this yields a complete PV representation of the interaction which we will call 'full 
PV representation in Fig. 1 and in the following. Such a decomposition can successfully 
describe the PV pion exchange on the Hartree-Fock level, i.e. the results calculated with 



the analytically known Hartree-Fock matrix elements ||2J] for the PV pion-exchange are 
reproduced. In the present formalism this can not be achieved by the other decompositions. 
On the other hand, however, e.g. the u exchange is no longer treated accurately in the 
full PV representation since PS admixtures arising from the Fock part of the u exchange, 
Eq. (|57|), are treated as PV. This bias will, however, turn out to be small. Anyway, all 
existing decomposition are unable to handle the PV pion exchange simultaneously with the 
remaining set of mesons on the Hartree-Fock level. As we will see later on, the one-pion 
exchange dominates the momentum dependence of the self-energy. Hence, it is reasonable 
to require that the PV pion-exchange is treated exactly on the Hartree-Fock level. This 
constraint is fulfilled adopting the full PV treatment. 



IV. RESULTS 



In the present section we study the impact of the various representations of the T-matrix 
(HI)' (13)' dHI) (|38D on the nucleon self-energy and on related quantities. 
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A. Momentum dependence of the self-energy 



On the level of the self-energy the effect of the different choices can be summarized as 

E(A;) = SPS(A;) - = SPS(A;) + ^ J -0^Mkq; qk) tr [GM{P - PV)] (39) 

with TF^ being the self-energy given in the pure PS representation. The different approaches 
for the self-energy S using for the T-matrix Eqs. (0), (|27|), (PH) or (^) are only varying 
in the choice of /r explained below. The strength of the suppression of the pseudo-scalar 
contributions is moderated by the amplitude which also determines the deviation (5S of 
the self-energy from the PS case. In principle this deviation shift is only apparent in the 
decomposition of S in the scalar and vector components, Eq. (|^), but vanishes when the 
complete matrix elements are built, i.e. 

< u{k)\5i:{k)\u{k) >= Q . (40) 

For the same effective mass M*{k-p) all approaches give the same total single particle po- 
tential ([l^) although they yield quite different scalar and vector self-energy components. 
However, the different approaches also yield different values for the effective mass M*{k-p) 
which leads to a different in- medium spinor basis \u > used in the self-consistent iteration 
procedure and the equivalence for the single particle potential gets distorted. Consequently, 
the different approaches lead to visuable changes also for those quantities which are built 
from total matrix elements, i.e. the single particle potential and the equation-of-state. If 
the value of M* is, however, kept fixed the equivalence on the level of matrix elements is 
exact and also numerically fulfilled with high accuracy. The cases discussed in the previous 
section modify only in Eq. (^) and can now be summarized as 

, PS 



R 



)^+^F5in-e) , 'PV choice' (41) 



MO) - f2{e) + h{e) , 'full PV 

Since earlier works on relativistic Brueckner theory which sticked to the projection method 
applied the scheme proposed by Horowitz and Serot (Eq. (|3ll) with P and P replaced 
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by PV and PV, respectively, our 'PV choice') we will discuss this method first. Within 
this scheme the structure of the self-energy, i.e. its density and momentum dependence 
has been investigated in detail in Ref. with the Bonn potentials as the bare NN- 
interaction. As the most prominent result we observed a strong momentum dependence 
of the scalar and time-like vector self-energy components around the Fermi momentum. 
These findings are in qualitative agreement with a previous work by Nuppenau et al. [Q. 
Fig. ^ shows the momentum dependence of the three self-energy components Eg, So, kpT^y 
at nuclear matter density p = 0.166 /m~^ obtained with the Bonn A interaction. The 
space-like Sv component is found to be relatively large inside the Fermi-sea and decreases 
with increasing momentum. Furthermore we compare to the corresponding reduced fields 



Ss, So where the S^-term is effectively included, Eq. (111). It is seen that the reduced 
fields are significantly smaller in magnitude at low momenta and generally show a less 
pronounced momentum dependence. The inclusion of the S^-term counterbalances the 
strong momentum dependence to some extent. The negative Sy-contribution effectively 



weakens the momentum dependence whereas in the opposite case ||6|JT9|| the momentum 
dependence will be enhanced. In the limit of a vanishing Sy-contribution reduced and 
projected components are equal. 

Fig.0 demonstrates the influence of the various mesons (using Bonn A). Taking only a 
and u exchange into account the result is quite similar to that obtained in Ref. i.e. 
the momentum dependence is flat inside the Fermi sea. Including the pion we are already 
very close to the full result. The strong momentum dependence of the present calculation 
originates to a large extent from the pion-exchange. Hence, the calculation still reflects the 
Hartree-Fock results (Fig.|l|) and the strong momentum dependence originates mainly from 
pseudo-scalar contributions of the pion. 

In Fig.^ the corresponding self-energies obtained for the various decompositions are com- 
pared. Adopting the full PV representation the space-like Sv contribution turns out to be 
much smaller than in the PS or the standard PV choice (see also Fig. 0). Therefore we show 
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the reduced self-energies Eg and So in which Ev is included for a better comparison. Both, 
the PS and the PV choice show a similar strong momentum dependence which again reflects 
the pseudo-scalar nature of the pion exchange. Consistent with the previous considerations 
this momentum dependence vanishes to a large extent when the pion contribution is sup- 
pressed by the full PV representation of the T-matrix. At high energies the different choices 
yield similar results since the influence of the pion decreases. The pure PS and the full PV 
representation can be regarded as the limiting cases which give the range of uncertainty in 
the determination of the self-energy. The full PV representation has thereby the big advan- 
tage that due to the weak momentum dependence the standard treatment of the Thompson 
equation, i.e. to approximate the effective mass by its constant value at the Fermi surface, 
can be safely applied as done in our and almost all other works on this subject. Furthermore 
this method ensures by construction that it treats the pion correctly at the Hartree-Fock 
level. Thus the full PV representation comes closest of all discussed representations to what 
we would call the optimal representation of the T-matrix. It is also worthwhile to mention 
that the results obtained within this representation agree well with recent calculations which 
include negative energy states and thus avoid the projection procedure |I9| . 

B. Single particle potential and the fit method 

As discussed above, the single particle potential, Eq. (p!3D , is in principle independent 
on the representation of the T-matrix which is again due to the on-shell equivalence of 
the corresponding matrix elements, Eq. (plf ). However, significantly different values of M* 
obtained in the different approaches lead to different results. Although this effect is reduced 
using the rescaled effective mass M*, the equation-of-state reacts sensitively on this M* 
dependence, as can be seen from table 1. A suppression of the PS contributions causes a 
larger effective mass and thus suppresses relativistic effects which originate from the mixing 
of small and large components of the spinors. A pure PS treatment leads to more binding and 
shifts the saturation point to higher densities. The corresponding equation-of-state is rather 
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stiff. By the empirical knowledge of the nuclear saturation properties the PS representation 
can be ruled out since the saturation density is much too high and the effective mass is 
much too small. The remaining two methods are in rough agreement with the empirical 
constraints, however, the densities are also slightly too high compared to the experimental 
Fermi momentum of about kp = 1.37/m~^. Here the larger value of the effective mass is 
a favor of the 'full PV representation since it is in better agreement with the constraints 
derived from the spin-orbit splitting in finite nuclei Compared to the 'PV choice', the 
'full PV representation leads to more binding and makes the equation-of-state softer, in 
agreement with the incompressibility derived from the isoscalar giant monopole resonance 
of about K ^ 230 MeV. One might assume that this is due to the fact that the 'full PV 
representation suppresses part of the repulsive ci;-exchange. This is however, not the case 
since the inaccuracy which arises in the PS i — > PV replacement procedure concerning the 
vector and tensor exchange amplitudes, Eqs. (0), has only a minor infiuence on the final 
results. We checked this point by treating the Hartree-Fock contribution of the one-pion- 
exchange separately in the PV representation while for the remaining part of the T-matrix 
the PS representation was retained. In terms of the self-energy this means to set fji = 
in (|55|) with the Hartree-Fock amplitude of the one-pion exchange. 



In view of the problems which arise in the determination of the self-energy we now 
briefly discuss a frequently used and much simpler approach first applied by Brockmann and 
Machleidt 0]. In this approach one tries to extract the self-energy components directly from 
the single particle potential, thus one avoids to take the explicit traces (|5|-^. Therefore there 
is no need to decompose the T-matrix into its Lorentz invariants. Actually, Brockmann and 
Machleidt determined only density dependent but not momentum dependent mean values for 
Ss, So by a fit to f/ given by Eq. (p!7|). This fit method works reasonably well as long as one 
restricts oneself to density dependent observables on the one-body level. The corresponding 
(reduced) effective masses are relatively close to our results at k-p obtained within the 'PV 
choice' [^. Thus it is understandable that the resulting nuclear matter saturation properties 
are similar in the two approaches of Refs. |^ and |^, see Tab.l. However, any information 
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on the magnitude of the space-hke Sv-contribution and the exphcit momentum dependence 
of the self-energy is completely lost in this approach. To overcome this drawback the Stony 



Brook group |T0| extracted momentum dependent fields from the single particle potential. 



They assumed a functional dependence of the (reduced) self-energies of the form 

at fixed density kp. The set of six parameters 7}s,o was then determined by a least 

square fit to U, i.e. by minimizing 



;^2 = y dJ.p 



( M* - 



H 2 



(43) 



However, such a procedure suffers from a large amount of arbitrariness since one tries to 
extract two independent functions out of one function. Consequently, the result is strongly 
infiuenced by the choice of the trial functions. It is clear that the class of trial functions 
given by Eq. ( ^2]) will generally not provide solutions of Eq. (p9|), although the insertion 
of the - a priori unknown - 'correct' amplitude should yield the correct results for the 
self-energies. 

To demonstrate this aspect we compare in Fig. | the fit procedure according Eqs. ( ^^5[ ) 
to the projection method using thereby both, the 'PV choice' and the 'full PV decomposi- 
tion. The fitted self-energies are determined from the single particle potential obtained with 
the 'PV choice'. It is seen that the fitted self-energies show a moderate momentum depen- 
dence which is in a qualitative agreement with the 'full PV representation, but not with 
the 'PV choice' to which the fit was performed. However, the asymptotic high momentum 
behavior of the fitted self-energies is completely different from the projected self-energies, 
independently which choice is used. Also the slope of the curves at low momenta seem to be 
distorted by the choice of the trial function. The results of Ref. [|10| show a similar behavior 
(which is due to the choice of the same functional ^-dependence (^)), however, the momen- 
tum dependence is a little more pronounced than our fitted results. For a fair comparison 
one has to be aware that in Ref. |jlO| and the present work different approximations to the 
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Thompson propagator were made. Actually, in the present work the Thompson equation 
is solved in the two-particle center-of-mass frame, whereas in Refs. PJl^ it is solved in the 
nuclear matter rest frame thus avoiding the projection techniques. The main difference is. 



however, that in Ref. ||T^ the effective mass entering into the Thompson equation and the 
effective spinor basis is allowed to be momentum dependent. This states an involved problem 
which affords a couple of additional approximations. If the momentum dependence of the 
self-energy is very pronounced one has to go beyond the present approximation scheme with 
M*{kp) in order to include this momentum dependence self-consistently. However, perform- 
ing the full PV decomposition (which we regard as the most reliable one) the momentum 
dependence is actually very weak (see Fig. 4). Thus the usage of a constant effective mass 



M*(|k| = kp, k-p) is well justified. In contrast to [|T0| where it was argued that the improved 



self-consistency suppresses the momentum dependence we find that momentum dependence 
is mainly governed by the representation of the T-matrix or - in the case of Ref. - by 
the choice of the trial functions. 

To illustrate this effect the single particle potential is considered in Fig. ^ It is seen that 
the different decompositions, 'PV choice' versus 'full PV representation, yield significantly 
different results for this quantity. The deviations are due to the different modifications of the 
effective interaction in the self-consistency scheme which yield also different effective masses. 
The 'full PV treatment lowers the potential by about 8 MeV compared to the standard 
treatment and thus leads to more binding in the equation-of-state. We also included the 
result of Ref. into this figure. Remarkably, this calculation yields the same result as the 
present 'PV choice' although the underlying scalar and vector self-energies are completely 
different. This somehow fortuitous agreement can be understood by the fact that the self- 
energies coincide around the Fermi momentum. In addition we show in Fig. ^ the single 
particle potential obtained as a result of the fit procedure, Eq. (^3]). U is reasonably well 
reproduced although the fitted and projected self-energies in Fig. 5 strongly differ. Hence, it 
is not possible to extract the self-energy components from the single-particle potential in a 
reliable way. This finding is also supported by a recent analysis where it was shown that the 
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fit method breaks down when apphed to isospin asymmetric nuclear matter [^]. Recently 
Miither, Ulrych and Toki |^ also determined momentum dependent scalar and vector self- 
energies components directly from the single particle potential (0). There M* was treated 
as an independent quantity which allows to generate more than one equation for f/(/c, k-p) 
to determine Eg and Sq. This approach, however, neglects that only the self-consistent M* 
has a physical meaning. With the method of ^ one also obtains a very weak momentum 
dependence of Eg and So similar to the 'full PV approach of the present work. 

Another experimentally accessible observable is the Schroedinger equivalent optical po- 
tential 0. Here the explicit momentum dependence of the self-energies provides an im- 
portant correction to the in first order linear energy dependence of the optical potential. 
Adopting the 'PV choice' we already found in Ref. [H a good agreement with the empirical 



values of Ref. for the real part up to energies around 800 MeV and an excellent agree- 
ment up to the pion threshold for the imaginary part. As can be seen from Fig. 7 where the 
real part of the optical potential is shown the agreement with the data is even improved at 
low energies when the 'full PV representation is used. This also indicates that the depth of 
the corresponding single particle potential is quite reasonable. At higher energies the two 
methods yield almost identical results. On the other hand, with the fields as obtained by the 
fit method and also predicted similar in [|l^ it is not possible to reproduce the high energy 
behavior of the optical potential. This will of course have implication when such fields, Eq. 
are applied to the transport description of heavy ion collisions. 



V. SUMMARY AND CONCLUSIONS 



In the present work we investigated the momentum dependence of the nuclear self-energy 
in the relativistic Brueckner approach. We applied the standard treatment which projects 
onto positive energy states and determines the self-energy components by a decomposition 
of the T-matrix into its Lorentz invariants. The T-matrix is represented by a set of five 
linearly independent covariants. Since the set of covariants is not uniquely determined in the 
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subspace of positive energies one has some freedom in the choice of the representation. It is 
therefore not possible to determine the scalar and vector parts of the relativistic nuclear self- 
energy in a unique way. This ambiguity originates from the fact that the on-shell scattering 
of positive energy states yields identical values for the pseudo-scalar and the pseudo-vector 
representation of the corresponding matrix elements and that they connect differently to 
the negative energy states, see also Refs. In the standard treatment of relativistic 

Brueckner theory one accounts for this fact by choosing a particular type of a pseudo- 
vector representation. To perform this 'PV choice', see Eq. (|3l|), one has to decompose 
already anti-symmetrized matrix elements into direct and exchange terms which is not free 
from ambiguities. Applying the Bonn potentials as the bare NN interaction we find that 
the conventional 'PV choice' leads to a pronounced momentum dependence of the nuclear 
self-energy. The spatial Sv-contribution of the self-energy is thereby found to be relatively 
large, in particular inside the Fermi sea, and counterbalances this momentum dependence 
to some extent on the mean field level. However, this strong momentum dependence makes 
the standard Brueckner approach questionable since the Thompson equation (or alternative 
reductions) are iterated using a self-consistent effective mass depending only on the Fermi 
momentum. 

The momentum dependence is found to be completely dominated by the pion exchange. 
It originates from the pseudo-scalar nature of the pion which is still remnant adopting the 
'PV choice' (|3T|) in the conventional manner. To eliminate this insufficiency we represented 
the T-matrix by a pure pseudo-vector decomposition and called this 'full PV representa- 
tion. The 'full PV representation accounts correctly for the desired pseudo-vector nature 
of the pion exchange and suppresses the momentum dependence of the self-energy almost 
completely. The two limiting cases, namely the full pseudo-scalar and the full pseudo-vector 
representation, set the range of uncertainty in the determination of the self-energy. How- 
ever, the 'full PV representation is more consistent with the usual approximation scheme of 
the Brueckner approach which assumes that the screening of the effective interaction in the 
medium introduces an additional density dependence, but is only weakly depending on the 
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momentum of the particles. We further investigated a frequently used method, namely to 
determine the scalar and vector self-energy components directly by a fit to the single particle 
potential. Although this method works reasonably well as long as one restricts to density 
dependent observables it leads to highly ambiguous results when applied to extract the full 
momentum dependence of the fields. Thus we conclude that the projection onto covariant 
amplitudes using thereby a complete pseudo-vector representation is up to now the most 
reliable way to determine the scalar and vector nucleon self-energy components as long as 
one works exclusively with positive energy states. 
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TABLE I. Binding energy per particle E, Fermi momentum kp, reduced effective mass M*, 
and compression modulus K for nuclear matter at saturation density employing the various repre- 
sentations of the T-matrix (41). As the bare NN interaction the Bonn A potential was used. 



fcp E M* K 

[fm-i] [MeV] [MeV] [MeV] 

PS 1.45 -17.70 455.6 335 

PV choice 1.41 -15.81 538.3 275 

fullPV 1.42 -16.59 648.8 245 

Ref. § 1.40 -15.59 564.3 290 
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FIGURE CAPTIONS 



Fig. 1: Pion contributions only to the nucleon self-energy in the Hartree-Fock approxima- 
tion. The self-energies determined in the pure pseudo-scalar (PS) or the pure pseudo- 
vector (PV, multiplied by a factor 10) representation are compared to results obtained 
within the 'PV choice'. In all calculations the nuclear matter density is chosen as 
p = 0.166 fm~^ and the pion coupling constant and form factor from the Bonn A are 
used. Solid lines represent the scalar, dashed lines the vector self-energy. 

Fig. 2: Momentum dependence of the nucleon self-energy in the relativistic Brueckner- 
Hartree-Fock approach using the 'PV choice' (see (^l|) . The nuclear matter density 
is chosen as p = 0.166fm~^ and the Bonn A potential is used. Solid lines represent 
the scalar, zero-vector and space-like vector components Sg, Sq and /cp^v; dashed 
lines represent the corresponding reduced fields into which the S^-term is effectively 
included (0). 

Fig. 3: Influence of the various meson exchange contributions on the nucleon self-energy 
(scalar part). The solid line corresponds to the full calculation (Bonn A), the dot- 
ted line to auo and the dashed line to crcuvr-exchange only. In all calculations the 'PV 
choice' is used. 

Fig. 4: Comparison of the reduced self-energy components, Eq. (|16|), obtained by the various 



decompositions of the T- matrix (see (^I]) . Solid lines correspond to the 'full PV 



representation, dashed lines to the 'PV choice' and long-dashed lines to a pure PS 
representation. The shaded area indicates the range of uncertainty spanned by the 
'full PV and the PS representation, Eq. (pT]) . The nuclear matter density is chosen 
as p = 0.166fm~^ and the Bonn A potential is used. 

Fig. 5: The reduced self-energy components obtained with the projection method using the 
'PV choice' (long-dashed, see (PD) are compared to the respective fields obtained by 
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a fit to single particle potential (^21) (dotted). In addition the results with the 'full 
PV representation (solid) and those of Ref. [jlO[ (dot-dashed) are shown. The nuclear 
matter density is chosen as p = 0.166fm^^ and the Bonn A potential is used. 

Fig. 6: Single particle potential in the relativistic Brueckner-Hartree-Fock approach with the 
Bonn A NN interaction at p = 0.166 fm~^. The results of the projection method with 
the 'full PV representation (solid hne, see (PD ) and using the 'PV choice' (dashed 



line) are shown. The dotted line represents the single particle potential obtained by 
the fit procedure (|42|) and (|43|) to the result of the 'PV choice'. For comparison also 
the calculation of Ref. [|10| is shown. 



Fig. 7: Schroedinger equivalent optical potential for Bonn A at p = 0.166 fm~^. The results 
of the projection method with the 'full PV representation (solid line, see (^1])) and 
using the 'PV choice' (dashed line) are shown. The dotted line represents the result 
obtained with the fit procedure (|4^) and (^31) to U. The empirical values (diamonds) 
are taken from pi 
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FIG. 1. 
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FIG. 2. 
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FIG. 3. 
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FIG. 5. 



33 



-40 



-60 



1 1 1 1 1 1 1 1 1 1 1 

fuUPV 


1 1 1 


PV choice 


* 


fit to PV choice 


0' 


□ Lee et al. 




// 
// 

- 

// 








1 1 1 1 1 1 1 1 1 1 1 




1 1 1 



0.0 0.5 1.0 1.5 

k [fm"^] 



FIG. 6. 



34 



100 - 



n ' 1 ' 1 

fuUPV 

PV choice 

fit to PV choice 

♦ Exp. 



50 - 







-50 - 




-100 



200 400 600 

k„-M [MeV] 



FIG. 7. 



35 



